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We compute 0(a 3 lna) relative corrections to the ground 
state hyperfine splitting of a QED two body bound state with 
different masses of constituents. The general result is then 
applied to muonium and positronium. In particular, a new 
value of the muon to electron mass ratio is derived from the 
muonium ground state hyperfine splitting. 
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The perturbative series for binding energies of a QED 
bound state are non-analytic in the fine structure con- 
stant and the expansion contains powers of a™ 1 In" 2 a, 
where ri\ and are some integer numbers. The ap- 
pearance of logarithms of a is best explained by the fact 
that different scales, such as the mass to, the momentum 
ma and the typical energy ma 2 control dynamics of the 
bound state. 

In a recent paper [Q we have explained how the non- 
relativistic Quantum Electrodynamics (NRQED) regu- 
larized dimensionally can be efficiently used to extract 
all In a corrections in a given order of the expansion in a 
and applied this technique to compute 0{a 3 In a) correc- 
tions to the decay rates of para- and orthopositronium. 
The purpose of this Letter is to apply that methods to 
the calculation of the C(a 3 lna) corrections to the hy- 
perfine splitting (hfs) of a general QED two body system 
in a ground state with an eye on the hfs in muonium and 
positronium. 

For both, muonium and positronium, there are good 
phenomenological reasons to consider 0{a 3 In a) contri- 
butions to the ground state hfs. The most precise mea- 
surement of this quantity in positronium gives [Q]: 

z£ s p = 203 389.10(74) MHz, (1) 

while the theoretical prediction which includes 

0(a 3 In 2 a) terms computed in ||, is: 

z/ t p h s = 203 392.01(46) MHz. (2) 

Obviously, the theoretical and experimental results dif- 
fer from each other by an uncomfortably large amount 
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(given the claimed accuracy of the two) , which indicates 
that further study of 0(a 3 ) corrections to the hfs of the 
positronium ground state is warranted. The complete 
calculation of 0{a 3 ) corrections is currently out of ques- 
tion because of tremendous technical difficulties; never- 
theless, the C(a 3 lna) corrections can be determined. 

Before considering the case of two equal masses, we de- 
cided to study a bound state of two particles with masses 
to and M . For to = m e and M = m M this corresponds to 
the bound state of the electron and the anti-muon called 
muonium. Our technique provides an excellent tool to ex- 
tract the C(o: 3 lna) corrections in this case keeping the 
full mass dependence and this result can be used in two 
ways. First, we will be able to check the correctness of 
our calculation against the known results for C(a 3 lna) 
corrections obtained in an expansion in m e /m^. Sec- 
ond, we will derive some new results from our formula; 
in particular, we will give the complete 0(a 3 m e /m fl In a) 
correction to the ground state hfs in muonium which is 
important for the extraction of the muon to electron mass 
ratio. 

The Letter is organized as follows. We first review 
our method of calculation (for more details we refer to 
fjj). We the n consider the unequal mass case (muo- 
nium) where virtual annihilation is not allowed. Later, 
we discuss phenomenological implications of our result 
for muonium and derive the new value of the muon to 
electron mass ratio. Finally, we describe how the calcu- 
lation should be modified in order to accommodate the 
positronium case and its phenomenological consequences. 

Our calculation is based on dimensionally regularized 
nonrelativistic QED with d = 3 — 2e being the number 
of spatial dimensions. In the NRQED framework two 
different contributions to the final result should be dis- 
tinguished. The first one is the hard contribution, which 
is sensitive to relativistic momenta only. This contribu- 
tion is not capable to produce any non-analytic depen- 
dence on a. The second contribution is the soft one. It is 
sensitive to nonrelativistic scales and for this reason can 
produce a non-analytic dependence on the fine structure 
constant. The main idea that permits a simple extrac- 
tion of the logarithmic terms is the following. In dimen- 
sionally regularized NRQED, the matrix elements of the 
nonrelativistic operators are the uniform functions of the 
fine structure constant. This implies that, when written 
in proper units, the dependence on a can be scaled out 
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of any matrix element. We refer to our recent paper 
for additional details on this approach; here we remind 
the reader that the relative momentum p, the relative 
coordinate r and the binding energy E scale as p — > jp, 
r — > 7~ 1 r, E — > j 2 jj,~ 1 E, with the scaling parameter 
7 = (>Za) 1/(1+2e) ||. Here fi = mM/(m + M) is the 
reduced mass of the bound state. We also assign the 
charge Z to the particle of mass M to distinguish recoil 
and radiative recoil contributions, as it is customary in 
bound state calculations. To illustrate how the scaling 
arguments help to compute the In a corrections, let us 
consider the matrix element of a nonrelativistic operator 
O that delivers 0(a 3 ) correction to the lowest order hfs 
kernel Veom (o" and £ are the spin operators of the two 
particles) , 
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We consider relative correction to the hfs: 
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The operator O is a function of r and p. Performing 
the rescaling of all the quantities on the left hand side of 
Eq.(EJ) according to the rules given above, we extract the 
dependence on a: 
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where n — 2, 3 is a power of a that explicitly enters the 
operator O and j,l,k are some integers. If the matrix 
element is finite, we can put e = and then the rela- 
tive correction to the hfs is a 3 times the a-independent 
function and hence no logarithms of a appear. There- 
fore, after the rescaling, the only place where In a can 
come from is the expansion of the factor ^ 3_n +J e in pow- 
ers of e; this implies that, in order to generate the In a 
corrections, the nonrelativistic matrix elements should 
diverge and only divergent pieces of the matrix elements 
are needed to determine the O(lna) corrections. Note 
also, that since the hard 0(a 3 ) contributions to the hfs 
are not needed, it is straightforward to keep the mass 
dependence exactly in our calculation. 

Let us now list all the contributions to the hfs in muo- 
nium and positronium relevant at C(a 3 lna). We begin 
with discussing the irreducible contributions, i.e. those 
that arise as average values of some local operators. The 
first of those is the operator that corresponds to a Taylor 
expansion of 0(a) hard scale contributions in powers of 
the relative momenta of the bound state constituents up 
to 0(p 2 ). For both radiative and annihilation corrections 
this contribution can be related to the divergences in real 
radiation [jjj . For the recoil contributions things are more 
complicated and the easiest way to extract the 0(ap 2 ) 
piece of the hard scattering amplitudes is to actually ex- 
pand the box diagrams to the required order. This results 
in the following contribution to the muonium hfs: 
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Other irreducible corrections at this order are pro- 
duced due to, loosely speaking, the magnetic moment 
renormalization of the 0(a 2 ) corrections to the hfs and 
in many cases the result can be simply obtained by gen- 
eralizing the calculation of Ref. [Q to the unequal mass 
case. We then derive: 
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for the contributions of the retardation and the "one- 
loop" operators, respectively (see [El for the nomencla- 
ture). Parameter £ = 1 distinguishes the contributions 
due to anomalous magnetic moment of the particle with 
the mass M. 

Two additional irreducible contributions originate 
from relativistic corrections to the single Coulomb or 
magnetic exchange when we account for the Pauli form 
factor in one of the vertices: 
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A M = - 
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The last set of irreducible corrections can be loosely 
described as the effect of the retardation on all the rele- 
vant operators that generate non-radiative corrections at 
lower orders. This includes: the third order retardation, 
the retardation in the one-loop operator, the irreducible 
exchange of two magnetic photons, the retardation in 
the graph with magnetic seagull operator on one line, 
the graph with two magnetic-Coulomb seagull vertices on 
both lines and, finally, the graph with the magnetic seag- 
ull vertex on each line. The resulting correction reads: 
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Several reducible contributions appear in the second 
order of time-independent perturbation theory. The first 
of them is generated by the so-called double seagull ef- 
fective potential. This contribution is very similar to the 
case of C(a 3 lno;) corrections to the positronium decay 
rate considered in [[TJ and it can be easily generalized on 
the unequal mass case: 
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The reducible retardation correction can also be de- 
rived as a simple generalization of the result in J5|]. The 
only difference is that the spin parts of the magnetic cur- 
rents also give non-zero contribution to the hfs. For un- 
equal masses, the result reads: 
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Also, the so-called ultrasoft contribution Q should be 
considered. We find: 
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For the last reducible contribution to the hfs, the 
corresponding nonrclativistic operator O is of the form 
BGVhi + Vh\GB, where G is the reduced Green function 
of the Coulomb Hamiltonian, B is the Breit Hamiltonian 
and Vhi is the hard radiative correction to B. The Breit 
Hamiltonian in d dimensions has been derived in [Q for 
M = m. For different masses, it reads: 
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where C is the Coulomb potential in d dimensions C (r) — 
-ZaT(d/2^1)/(TT d ^- 1 r d - 2 ). 

The potential Vhi is the sum of four contributions ||: 
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where the first one arises from the Coulomb photon ex- 
change with one of the vertices being either the one-loop 
slope of the Dirac form factor or the Pauli form factor, 
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the second one is due to the one-loop anomalous magnetic 
moments, 

the third one comes from the hard one-loop box dia- 
grams, 
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and the last one accounts for the one-loop vacuum polar- 
ization (hadronic vacuum polarization is not included): 
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The structure of Viii is very similar to the structure 
of the Breit Hamiltonian and so the corresponding cal- 
culation goes along the lines of M. On this way one 



recognizes that both D wave and S wave contributions 
should be considered. 

Since the D wave part of the new perturbation (lq) 
differs from that of the original Breit perturbation (13|) 
only by the overall factor, we can read off the D wave 
contribution to the 0(a 3 In a) correction to the hfs from 
the corresponding 0(a 2 ) correction in positronium (see 
I): 
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To find the contribution of the intermediate S states, 
we first project both B and Vhi on to the S wave and 
then proceed along the lines described in Q). Since all 
the steps in this calculation have their counterparts in the 
calculation described in detail in Q and since the inter- 
mediate formulas for the different mass case are lengthy, 
we refrain from presenting them here. 

Summing up all the relevant contributions, we obtain 
the final result for the 0(a 3 In a) correction to the hfs of 
the ground state in the unequal mass case (we put £ = 1 
below) : 
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Eq.(gOj) is one of the principal results of this Letter. 
If we identify m with the electron mass and M with the 
muon mass, Eq.(^o|) provides the result for the 0(a 3 ln a) 
correction to muonium hfs. Muonium has been studied 
extensively over the years and much is known about this 
system. In particular, there are certain limits of Eq.([20|) 
that can be checked against known results. To this end, 
it is instructive to expand Eq.(pO|) in powers of m/M up 
to the first non-trivial order: 
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The last equation shows that our result, Eq.(^o|), 
correctly reproduces the \n 2 (Za) terms as well 

as the Z 2 \n(Za) single logarithmic term [7@j and the 
{m e /m^)Z 3 \n(Za) ln(m^/m e ) term |9|,|irj|| which are all 
available in the literature. 

We now proceed to the discussion of what this result 
implies for the phenomenology of the muonium hfs. We 
first note that in this case Eq.(21) can be used since 
higher powers in the expansion in m^jm^ have a negli- 
gible impact. Since the ln(m^/m e ) enhanced part of the 
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0{m & lm tL ot'\w.a) corrections has been properly taken 
into account in a recent compilation of all theoretical re- 
sults for muonium hfs we disregard it here. Setting 
Z = 1, we then obtain the 0(m e /m tl a i In a) hfs shift: 
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Numerically, it evaluates to 0.013 kHz; this should 
be compared with a similar contribution <5V Mu = 
—0.265(64) kHz, originating from the incomplete calcula- 
tion in , that has been accounted for in the theoretical 
value for muonium hfs in @. 

The difference between the two numbers has significant 
impact on the electron to muon mass ratio determination. 
It is easy to see, that it amounts to the relative shift of 
6.2 x 10~ 8 in this ratio (compare with the quoted relative 
theoretical uncertainty 2.7 x 10~ 8 |12|) and, if we use the 
central value from |l2|, Eq.(161), the new result reads: 



— ^ = 206.768 2784(30)(23). 
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Here the first error is related to the error in the theo- 
retical prediction for the muonium hfs and the second is 
the experimental one. The theoretical error in the hfs was 
estimated following Jl^] ; the only difference is that we es- 
timated the uncalculated non-logarithmic Ofae/m^a 3 ) 
recoil and radiative- recoil corrections as half of their In a 
enhanced counterparts in Eq. (^l|) . The uncertainty in 
the muonium hfs due to uncalculated higher order correc- 
tions we have obtained in this way is 0.07 kHz, compared 
to 0.12 kHz in @. 

For positronium, the calculation goes essentially un- 
changed, although two facts have to be noticed. First, 
since there is an annihilation contribution to the leading 
order hfs in positronium, the relative weight of different 
corrections changes. Second, one has to take into ac- 
count additional annihilation contributions to the Brcit 
and Vhi operators. These annihilation operators read (S 
is the spin of positronium): 



positronium ground state hfs, so that the theoretical pre- 
diction becomes: 



f t p h s = 203 391.69(16) MHz. 
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The theoretical prediction moves closer to the experimen- 
tal result in Eq.([T]), however the difference is still signif- 
icant. Since the value of the new 0(ma 7 In a) contribu- 
tion turns out to be roughly one third of the Oima 7 In 2 a) 
one, the series look reasonably convergent. For this rea- 
son we estimate the nonlogarithmic 0(a 3 ) contribution 
to the positronium hfs as being one half of the logarith- 
mic one. This is the origin of the uncertainty estimate in 
Eq.©. 

In conclusion, we have computed C(a 3 lna) correc- 
tions to the hyperfine splitting of the general QED bound 
state keeping the full mass dependence. We then applied 
this result to the hfs of the muonium and positronium. 
The new value for the muon to electron mass ratio is 
extracted from the ground state hyperfine splitting in 
muonium. As for the positronium ground state hfs, the 
computed correction slightly reduces the discrepancy be- 
tween theory and experiment. However, it is hard to 
imagine that higher order corrections can further signifi- 
cantly shift the theoretical value. In this circumstances, 
one should perhaps start taking the discrepancy between 
the theory and experiment in the positronium ground 
state hfs seriously. 
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and they should be added to B and Vhi, respectively. 
Finally, for obvious reasons, one should disregard the 
0(1/M 2 ) contribution in Eq. (fL8|) . Proceeding along the 
lines described for the unequal mass case, we derive the 
result for the 0(ma 7 ) hfs in positronium: 
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Numerically, the new 0(ma In a) term gives an addi- 
tional shift of —0.32 MHz to the theoretical value of the 
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